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SUMMARY 

A laser speckle, differential strain measurement system has been built 
and tested for future applications In hostile environments. One-dimensional 
electronic correlation of speckle pattern movement allows a quasi-real time 
measure of strain. The system has been used successfully to measure uniaxial 
strain reaching Into plastic deformation of a test specimen, at temperatures 
ranging to 450 °C. A resolution of 16 microstrain Is given by the photodiode 
array sensor pitch and the specimen to sensor separation. The strain meas- 
urement error Is estimated to be ±18 microstrain ±0.3 percent of the strain 
reading. The upper temperature limit of the gauge Is determined by air density 
perturbations causing decorrelation of the reference and shifted speckle pat- 
terns, and may be Improved by limiting convective flow In the Immediate 
vicinity of the test specimen. 


INTRODUCTION 

This paper covers the Phase I results of a multiphase effort to develop an 
optical strain gauge capable of mapping the principal axes of strain on the 
surface of a hot specimen. Durability enhancement programs for propulsion 
system components prompted the need for such Instrumentation. The objective of 
this first phase has been to provide a noncontact, single axis, differential 
strain gauge for experimental purposes. Sverdrup Technology, Inc. Is develop- 
ing this optical strain measurement system for the Lewis Research Center. 

The design of the optical strain gauge Is based on the speckle shift 
method of Ichlrou Yamaguchl^ at The Institute of Physical and Chemical 
Research, Japan. Laser speckle Is a phase effect caused by the diffuse reflec- 
tion of spatially coherent light off of a rough surface. Interference of the 
light reflected off of a specimen results In a light and dark Intensity distri- 
bution at the sensor. Yamaguchl's noncontact method of measuring differential 
surface strain makes use of speckle field displacements detected by cross- 
correlating electrical signals from a linear photodiode array. A limiting 
geometry makes speckle displacement due to stress on a test specimen directly 
proportional to the associated strain. Use of this method eliminates problems 
usually caused by subjecting conventional measurement devices to hostile envi- 
ronments such as extreme temperatures, and high electromagnetic fields; the 
sensor need not be In the same environment as the specimen. In addition, 
Yamaguchl's method has the advantages over other optical strain measurement 
techniques of not only being noncontacting, but needing no surface preparation, 
automatically correcting for strain error due to rigid body motion, and giving 
quasi-real time results (the value of strain Is calculated within seconds of 
taking the data). 
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THEORY 


The utilization of laser speckle to determine differential strain Is 
achieved by applying the Fresnel-KI rchhof f diffraction Integral describing the 
propagation of light In free space. This Integral Is used to derive the rela- 
tionship between surface deformation of a specimen subject to loading and 
speckle displacement In the diffraction field (refs. 1 and 2). The speckle 
displacement contains terms of translation, rotation, and strain, from which 
the principle strain term must be extracted. This extraction relies on the 
use of two symmetrically Incident laser beams reflected sequentially onto a 
linear photodiode array oriented parallel to the specimen surface. After the 
reference (before-strain) and shifted (after-strain) speckle patterns from 
each beam are recorded and correlated, the difference In speckle shift between 
them Is taken leaving only the component of shift due to surface strain. 

The theory Interpreting the speckle displacement on the observation plane 
Is based on the assumption that the Illuminated region of the specimen Is homo- 
geneous; this region acts as a random diffraction grating. For speckle shift 
to accurately represent strain then requires that the change In the surface 
roughness spacing be uniform. This proved to be the case In the experimenta- 
tion within the elastic limit of the test specimen, as observed In the stabil- 
ity of the speckle patterns and the agreement of the results with known 
quantities. 

Referring to the simplified coordinate system shown In figure 1, a laser 
beam Is Incident on the diffuse surface of a test specimen at an angle e s 
from the normal to the surface. The observation plane Is the X,Y plane, and 
the object plane Is the lowercase x,y plane. The center of the Incident laser 
spot Is at the origin of the x,y plane. When the object Is stressed and 
deformed, the object point Is displaced bylftx.y). The ra<J1us of curvature of 
the wavefront Is denoted L$. The speckle displacement A Is broken down 
Into components, given by 

A X 3 a x ' L 0 [ c xx s1n( V ” fi y [ C0S ( e s ) + 1 ( 1 -a) 

A y = a y - l Q j’c xy s1n(e $ ) - Q x [cos(e $ ) + 1) - ft z s1n(e $ )j (l.b) 

where ~t, Q, and [T] are the translation vector, the rotation vector, and stain 
tensor, respectively. As long as the displacement Ay Is much smaller than 
the sum of the sensor height and the typical speckle size on the sensor Ax 
can be determined by cross-correlating the sensor signals before and after sur- 
face deformation. The correlation Is performed between the reference and the 
shifted speckle patterns from the left Incident beam, and again between the 
reference and the shifted patterns from the right Incident beam. 

Initial cancellation of rigid body motion terms requires the Incident radius 
of curvature L 5 to be much greater than the sensor distance Lq, and Lq to 
be much greater than the spot diameter of the Incident beam. Equations (1) 
show the remaining rigid body motion components In Ax, which can be cancelled 
out by taking the difference between speckle displacements from the two beams 
at +e s and -e$. From equation (l.a) then 


2 



AAx = -2LQe xx s1n(0) (where 6 = 1 65 1 ) 


( 2 ) 


or, 


. - aA x 

e xx ” 2L Q s1n(e) 


(3) 


If the criterion Lq/L$ « 1 Is not met, then the out-of-plane error term 


error term = 


-a z cos( 6 ) 


(4) 


Is added to the strain formula (eq. (3)). 


MEASUREMENT SYSTEM DESIGN AND SETUP 

The system Is designed to run predominantly under computer control, with the 
operator selecting procedures from menus on the CRT. A schematic of the opti- 
cal setup Is depicted In figure 2. The laser and the beamsteering components 
are mounted on a hard-coupled optical table (this non-lnterferometrlc tech- 
nique does not require vibration Isolation). Two variable beamsplitters are 
used to create three beam paths, each beam switched by an acousto-optic modu- 
lator. A line scan camera records the speckle patterns generated by the left 
and right beams at 6 $ = 45°, and an area array camera allows the test section 
and gauge position to be monitored remotely. The cameras, modulators, and the 
motorized mirror mounts are remotely controlled. The testing machine and 
Induction heater are mounted on a hydraulic lift table to provide positioning 
relative to the optics. An error analysis of this system shows an overall 
strain measurement error estimate of ±18 microstrain ±0.3 percent of the strain 
reading. 

The control and data paths are shown In figure 3. The system controller 
Initiates procedures through three main components, these being the data 
acqulsltlon/control unit (DAQ), and two pulse generators. The data acquisition/ 
control unit outputs control signals to the custom circuitry used for exposing 
and reading the linear and area array cameras; these circuits synchronize the 
pulse generators with the array cameras, triggering the pulse generators to 
output TTL exposure pulses to the acousto-optic modulator RF generators. The 
triggered outputs of the pulse generators switching on and off the beam- 
deflecting modulators act as high speed programmable shutters. The DAQ Is 
also used to read the voltages from thermocouples, resistance strain gauges on 
the specimen, and a load cell. Higher temperatures are also measured by the 
Infrared pyrometer which, In addition, provides a temperature control signal 
to the Induction heater. This allows precise control of the specimen tempera- 
ture at the gauge position. 


RESULTS AND RECOMMENDATIONS 

A series of test results shown In the following paragraphs give a good 
Indication of the overall system response. A 750 diode cross section of each 
speckle pattern Is used In the correlations In order to minimize the execution 
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time of the measurement. A typical pair of speckle patterns (reference and 
shifted) and their cross-correlation function Is shown In figures 4(a) and (b). 
The shift between the speckle patterns was Induced by a 13.8 MPa (2 ksl) load. 
The shift A x was calculated to be 7 diodes (105 ym) , corresponding to 
65 microstrain on the specimen after rigid body shift was cancelled out. The 
correlation curve was not normalized. In order to decrease the execution time 
Involved. 

Some plots of typical data runs are shown In figures 5 to 8. These figure 
show data for Isothermal runs at room temperature and 450 °C, as well as a run 
measuring thermal strain. All optical data shown here were taken by summing 
Incremental shifts between each point. A load offset was applied to the speci- 
men before each run; therefore, the stress axes do not represent an absolute 
value of stress. The test specimens are made from Inconel 600 stock, an alloy 
used for high temperature applications. 


Results 

Tests conducted during this program verified the stability of speckle pat- 
terns generated by an oxidized specimen subject to maximum strains over 
0.1 percent. The stability of the surface structure of the specimen was not 
observed to be affected by high temperature once oxidation occurred. The tests 
were conducted In an open air atmosphere. 

Figure 5 shows the optical gauge data plotted alongside the data from a 
resistance strain gauge taken at the same time. The curves can be seen to 
enter the nonlinear region of strain, and their slopes agree to within 6 percent. 
An offset of 6.9 to 13 MPa (1.0 to 2.0 ksl) occurs at the beginning of the run, 
where no strain Is measured optically. This offset (Initially large slope) Is 
typical of the optical measurements at the start of the load cycle, and Is 
believed to be caused by movement of the specimen at small angle deflections 
of the eccentric crank connected to the load arm of the testing machine. 

In figures 6(a) and (b), hysteresis due to plastic strain Is observed In 
both the optical and resistance gauge data as the load Is first Increased to 
262 MPa (38 ksl) and then backed down to the Initial preset value. Both sets 
of data were generated on the same run. Figure 6(a) shows the plots of the 
optical gauge data, and figure 6(b) shows the corresponding resistance gauge 
data. Agreement between the two gauges Is very good, with the final offset 
strain values of the optical and resistance gauges within 2 microstrain of each 
other. This Is well within the resolution of both gauges. Curve fitting over 
the linear portions of the data show differences of 8 percent on the load up 
cycle and 3 percent on the load down cycle. 

The optical data curve Is discontinuous at the beginning of the load down 
cycle, due to the out-of-plane displacement term a z . This type of break In 
the curve Is often observed as the load Is relaxed and the specimen returns to 
Its neutral, unloaded position along the z axis. In figure 6(a), the a z 
offset on the load down cycle Is nearly equal to the sum of the a z offsets 
In the opposite direction during the load up cycle. The remaining a z dis- 
placement Is spread smoothly along the curve as a function of load. It Is 
observed that when the plots from each gauge are superposed and offset to 
align the load down curves, the slope signatures are the same. 
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Figure 7 shows stress-strain data at high temperature; the temperature 
setpoint was 450 °C. The straight line drawn on the plot Is the modulus of 
elasticity E(T = 450°) given by an Inconel 600 data sheet. The slope calcu- 
lated from the optical gauge data Is 10 percent higher than the handbook value. 
The optical data show discontinuities below 67 MPa (10 k s 1 ) due to specimen 
shifting In the z direction (a z ), but smooths out for a constant slope up to 
200 MPa (29 ksl). At high loads, the term a z has been observed to Increase 
proportionally with load, which accounts for the slightly larger slope of the 
optical data. 

A test run measuring thermal strain Is shown In figure 8. As with the Iso- 
thermal runs, poor correlation occurred at temperatures above 450°, and the run 
was terminated. A linear fit of the data gives the measured coefficient of 
expansion to be 14.3 ym/m’C. This Is within 3 percent of the handbook value. 
Some random error (±14.8 pm/m) occurred due to the uncertainty In the temper- 
ature setpoint value. 

A least-squares regression over the linear portion of the strain data taken 
for each run provides a good measure of the accuracy of the data. The calcu- 
lated modulus of elasticity (E) or coefficient of thermal expansion shows error 
trends In the system (e.g., the error due to a z Increasing with load). The 
linear-correlation coefficient Indicates the degree of confidence In E, or how 
well strain Is shown to be a function of stress. The plots shown here have 
linear-correlation coefficients ranging from 0.998 to 1.000. 


Phase I Conclusion and Phase II Design 

The experimentation conducted during this project has resulted In a number 
of useful Insights regarding the Implementation of Yamaguchl's theory. Over- 
all, this measurement system works very well for measuring strain at low to 
midrange temperatures (up to 450 °C) using an open air specimen configuration. 
Increasing the temperature range depends on how well air density perturbations 
In the vicinity of the specimen can be controlled. These perturbations cause 
time-varying speckle Intensity fluctuations which lead to a relative shift 
between the left and right correlation functions. The temperature range could 
be Increased by placing baffles In the vicinity of the heated specimen, or by 
enclosing the specimen, which would block convective air flow and decrease 
transverse thermal gradients between the specimen and the sensor. 

The most practical result of the testing, though, Is that the effectiveness 
of Yamaguchl's speckle shift technique as a strain gauge Is limited by the 
extent of the system's Immunity to rigid body motion of the test specimen. 

Thus, the limiting factors In Implementing the test system are terms of rigid 
body motion. An Increase In the radius of curvature of the laser beam wave- 
front, as well as careful alignment of the specimen In the load train, can 
eliminate the systematic error observed In these tests and ensure the reliabil- 
ity of the data with minimum equipment specifications. The present system 
contributes an out-of-plane error term (eq. (4)) of 0.13 a z strain to the 
measurement. The a z sensitivity can be decreased sufficiently by Increasing 
the radius of wavefront curvature L$. For the case of a test situation where 
translations of the test components are extreme and unavoidable, the a z 
translation term could be measured for each pair of speckle patterns recorded 
and calculated out of the speckle shift equation. 
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The other effect of rigid body motion Is decorrelation of the reference and 
shifted speckle patterns. This Is caused primarily by a prohibitively large 
component of shift In the vertical direction (Ay). As Ay Increases, the 
amount of common Information between the reference and shifted speckle patterns 
as recorded by the photodiode array becomes less and less. If Ay Is a slowly 
varying function with respect to load, then decorrelation can be avoided by 
Incrementing the load In small steps, and summing the strains as the reference 
patterns are updated. The sensitivity to Ay can also be decreased by placing 
a cylindrical lens between the specimen and the linear photodiode array. 

The second phase of the Optical Strain Measurement System Development Is In 
progress. A schematic of the proposed Phase II system Is shown In figure 9. 

The Phase II system design allows the plane of Incidence to rotate about the z 
axis. Since the Intersection of the plane of Incidence with the specimen sur- 
face defines the sensitive axis of the gauge, rotation of the plane of Incid- 
ence enables multiple components of strain to be measured and the principal 
strain axes to be calculated. A single beam Is switched between the two Incid- 
ent paths by a Pockels cell and polarizing beamsplitter cube combination. 
Rotational Invariance between the fixed and rotating assemblies Is achieved by 
a pair of 1/4 wave plates. The radius of curvature L 5 of the Incident beam 
Is Increased by two orders of magnitude by a lens system which positions the 
beam waist at the specimen surface. 
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FIGURE 2. - OPTICAL SETUP AND TEST EQUIPMENT. 
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FIGURE 7. - DATA FROM THE OPTICAL GAUGE AT R50 degrees C. 



FIGURE 8. - OPTICAL MEASUREMENTS OF THERMAL EXPANSION. THE STRAIGHT LINE INDICATES 
THE HANDBOOK COEFFICIENT OF THERMAL EXPANSION. 
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